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We derive the Gauss-Codazzi equation in the holonomy and plane-angle representations and we use 
the result to write a Gauss-Codazzi equation for a discrete (2-|-l)-dimensional manifold, triangulated 
by isosceles tetrahedra. This allows us to write operators acting on spin network states in (2-1-1)- 
dimensional loop quantum gravity, representing the 3-dimensional intrinsic, 2-dimensional intrinsic, 
and 2-dimensional extrinsic curvatures. 


I. INTRODUCTION 


Three distinct notions of curvature are used in general relativity: the intrinsic curvature of the spacetime manifold, 
M, the intrinsic curvature of the spacial hypersurface E embedded in M which is utilised in the canonical framework, 
and the extrinsic curvature of E [1]. The three are related by the Gauss-Codazzi equation. On a discrete geometry, 
the definition of extrinsic curvature is not entirely clear [2]. Even less so in loop quantum gravity (LQG), where 
the phase-space variables are derived from the first order formalism [3], which generically does not yield a direct 
interpretation as spacetime geometry. But a definition of these curvatures is important in LQG, because we expect a 
discrete (twisted) geometry to emerge from the theory in an appropriate semi-classical limit [4]. 

In n-dimensional discrete geometry, the manifold is formed by n-simplices [2, 5]. The intrinsics curvature sits on 
the (n — 2)-simplices, called hinges, and can be defined by the angle of rotation: a vector parallel-transported around 
the hinge gets rotated by this angle. The rotation is in the hyperplane dual to the hinge. 

In the canonical formulation of general relativity it is convenient to use the ADM formalism [6] or a generalisation 
[7]. The phase space variables are dehned on an n — 1 surface: the “initial time”, or, more generally, “boundary” 
surface E. At the core of the ADM formalism is the Gauss-Godazzi equation, relating the intrinsic curvature of the 
4D spacetime with the intrinsic and extrinsic curvatures of E. In the context of a discrete geometry, E is defined 
as an (n — I)-dimensional simplicial manifold formed by (n — l)-simplices. To develop a formalism analogous to the 
ADM one, we need an equation relating the intrinsic and extrinsic curvatures of E and M, a ’discrete’ version of the 
Gauss-Codazzi equation. 

In this paper, we explore a definition of extrinsic curvature for the discrete geometry which allows us to write an 
equation relating the extrinsic and intrinsics curvatures, which we refer as the Gauss-Codazzi equation for discrete 
geometry. 

To get some insight into this problem, we hrst review the standard Gauss-Codazzi equation, both in second and 
first order formalism. In Section II we derive the holonomy (matrix) representation and plane-angle representation 
of the Gauss-Codazzi equation. We consider the discrete version of the Gauss-Codazzi equation in Section III, where 
we define all the curvatures. In Section IV, we move to the loop quantum gravity picture and study the operators 
corresponding to these curvatures into operators. 


II. GAUSS-CODAZZI EQUATION 
A. Standard Gauss-Codazzi equation 


The Gauss-Codazzi equation relates the Riemann intrinsic curvatures of a manifold and its submanifold with 
the extrinsic curvature of the submanifold. We will briefly review the continuous (2-|-l)-dimensional Gauss-Codazzi 
equation in this section, but the formula will be valid in [n + l)-dimension. 
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1. Second order formulation 


Consider a 3-dimensional manifold M and a 2-dimensional surface S embedded in M. In the second order formu¬ 
lation of general relativity [8], the Riemann intrinsic curvature of M is a ( „ ) tensor: 


= ^R^yjsda A dx^ ® dx^ A dx'", a, = 0,1, 2 

which can be thought as a map that rotates a vector v G TpM parallel-transported along an infinitesimal square loop 
defined by unit vectors Next, we define the projection tensor as: 

Qui' = - (n,n)n^n^, 

with Up is the normal to hypersurface S C M. {n,n) = ± 1 , depends on the signature of the metric. The projected 
Riemann curvature of on S is defined as: 

^R\j. = A dx^ ® dx^ A dx" = q^ Qp ^R’^'v'p'^a A dx^ ® dx^ A dx''. 

Taking the projected part from the full part we have relation as follow: 

+ S, (decomposition formula), (1) 

S is the ’residual’ part of ^R. The projected part ^R\^ can be written as: 

^R\^ = ^R+ [K,K] (Gauss equation), (2) 

which can be written in terms of components: 

"K.p = - KKpp), ( 3 ) 

^R and K are, respectively, the 2-dimensional Riemannian curvature and extrinsic curvature of hypersurface S. The 
residual part S can be written as: 

(Codazzi equation), 

with is the covariant derivative on the slice S. 

Writing (1) and (2) together, we obtain: 


^R= ^R+[K,K]+S. 




2. First order formulation 

The gravitational field is a gauge field and can be written in a form closer to Yang-Mills theory. This way of 
representing gravity is known as first order formulation of general relativity [3, 8]. Let manifold M be a spacetime. 
Let e be a local trivialization, a diffeomorphism map between the trivial vector bundle M x with the tangent 
bundle over M: TM = Up {p x TpM}^ and A be the connection on M x The 3-dimensional intrinsic curvature of 
the connection is the curvature 2-form: 

® dx^ A dx'', fi, v,I,J = 0,1,2. 

which comes from the exterior covariant derivative of the connection: 

F = doA. (4) 

{dp,dx^} and {^ 7 ,^^} are local coordinate basis on M and respectively. In terms of components: 

RIuj = 9fj,Alj — dtjApj + Apj^A^j — Alj^Apj 


(5) 
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We use nice coordinates such that the time coordinate in fibre is mapped by e to the time coordinate do in the 
basespace M, i.e., the map e is fixed into: 

e (^0) = e° = = dx°. 

For the foliation, we use the time gauge, where the normal of the hypersurface S is taken to be the time direction: 

n = do, (n) = ^o- 

Then the (2 + 1) split is simply carried by spliting the indices as / = 0, a, and /i = 0, i, with 0 as the temporal part: 

= { — FaOiO + Fooia + ^oiOO ~ FoiOa) A 0 dx^ A + {FaijO ~ Fmja) A 0 dx^ A dx^ + (6) 

+ (Faiob - FaOib) A (g) dx'- A dx° + Faijh^^" A (g) dx* A dx^ 

The projected part of the curvature 2-form on S is: 

= FaijbC A (g) dx* A dx\ 

and the residual part is: 

s = (-FaOio + Fooia + Fjiioo “ ^oiOa) ^ A ® dx* A dx° -f {FaijO - Fo,ja) r A (g) dx* A dx^ 

+ {Faiob - FaOib) c A (g) dx* A dx°. 

Therefore, the decomposition formula is clearly: 

^F= ( 7 ) 

Next, we take only the projected part: 

^F\^ = A (g) dx* A dx-’, 

which can be written in terms of components as follow: 

Fa^Jb = "^Faijb + (Gauss equation). (8) 

The closed part of ^F\^ is clearly the 3D intrinsic curvature of connection in E and the rest is the extrinsic curvature 
part: 

^F = ^Fa^jbC A G dx* A dx^ = {d^Aa^b - d,Aa,b + Aa.^A^^^ - AajaAl) r A e"' G dx* A dx^ 

K = Kfia. ® dx* = A%i^ dx\ 

The residual part S, in this special coordinates and special gauge fixing satisfies: 

Foija — FjKia F^Kja (Codazzi equation), 

and the other temporal components are zero. It must be kept in mind that this is the Gauss-Codazzi equation for a 
special local coordinate and special choice of gauge fixing, for a general case, they are not this simple. 

To conclude, we have the decomposition and Gauss-Godazzi equation for a fibre bundle of gravity: 

=‘^F+ [^K, "^K] +S. 

B. Gauss-Codazzi equation in holonomy representation 

1. Holonomy around a loop 

In this section, we will write the Gauss-Godazzi equation in terms of holonomy. The holonomy is defined by the 
parallel transport of any section of a bundle, say, sq, so that it satisfies the equation as follow: 

Dv (so) = dso + A (sq) = 0. 


( 9 ) 



Solving (9) using recursive method [3, 7, 8 ], we obtain the solution; 

s{t) = U (7 {t),D)so, 
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U (7 (t), D) is holonomy of connection D along path 7 (t): 


U{-f{t),D) = U^ = Pex-pjA, (10) 

with P is the path-ordered operator (See [3, 7] for the details of the derivation). 

Consider a square loop 7 embedded in S which encloses a 2-dimensional area. The holonomy around the square 
loop can be written as a product of four holonomies, since holonomy is piecewise-linear: 

^ 

Taylor expanding the holonomy in (10) up to the second order [3, 8 , 9], we obtain: 

U-f = \ -\— —F^v + O , 

with is an infinitesimal square area inside loop 7 . The formula: 

= U^-l-0^ ( 11 ) 

will be used to write the Gauss-Codazzi equation in terms of holonomies. 


2. First order formulation in holonomy representation 
Contracting ( 6 ) by an infinitesimal area 0 ^“^, we obtain: 


^iO 


Fli_tuJ = ^UlJ — SlJ — = — { — FaOiO + Tboia + FaiOO — FoiOa) S + {FaijO — Fenja) A -I- 

+ ^ {FarOb - Fa0^b) T A A C', 

where we have used the result in (11), relating the holonomy with the curvature of the connection. Taking the 

projected part ( 8 ) contracted with a®-^, and using the relation between holonomy with the curvature 2 -form, we can 

write the projected part in terms of 2 -dimensional holonomy and the contracted extrinsic curvature as follow: 


-YFa^Jb = -Y ^Fa^Jb + -Y {Ka^KJb " K^jK^b) = ^Uab - Sab - + — {Ka^K,b - KajK.b) ■ 

Finally, collecting all these result together, we obtain: 

-Y- ^ = T A ^ {Ka.K.b - KajK^b) C A ^ 


+ ^ (( —-FaOiO + Tboia + FaiOO — FoiOa) A ^° -|- {FaiOb — FaOib) A 

+ ^{Fa^JO-Fo^,a)£.’^ Fe 




or simply: 


+[K,K] +S + - ^I), (12) 

'-,,-^ 

T 

where we have written the Gauss-Codazzi and the decomposition formula together, noting the residual terms as S. 
Remember that U is the holonomy of the connection, which describes the curvature of the connection of fibre F, not 
the curvature of the basespace M. 
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3. Second order formulation in holonomy representation 
In the same manner as above, for a tangent bundle TM, we obtain the relation 

-Y Adx^ = + Y ^Rkijidx'" Adx’- + — {KkiKji - K^jK^) dx^ A dx^ 


/\dx^ {^Hki—Ski—0‘^)dx^Adx’' [K,K]f^^dx^ Adx'- 

H“ ( ^jOiO RoOij RjiOO RoiOj^ “1“ {RjikO ^Oikj^^ dx^ A dx 
+ ^ {RjiOk - RjOik) dx^ A dx’^ 

or simply: 

= ^H+[K,K] + S+{^I- 2/), 


5 , 


(13) 


with H is the holonomy around any loop 7 embedded in S C M, coming from the Riemann tensor R. 


C. Gauss-Codazzi equation in plane-angle representation 

Rotations can be represented in two ways, using the holonomy (matrix) representation, or using the plane-angle 
representation. For a continuous theory, the holonomy representation provides a simpler way to do calculation 
concerning rotation. In a discrete theory, where we would like to get rid of coordinates, the plane-angle representation 
is a natural way to represent rotations and curvatures. The plane-angle representation is only a representation of the 
rotation group using its Lie algebra for the plane of rotation and one (real) parameter group times the norm of the 
algebra for the angle of rotation. There exist a bijective map sending the holonomy to the plane-angle representation, 
known as the exponential map [10]. 

In this section, we rewrite the Gauss-Codazzi equation using the plane-angle representation. We only do the 
calculation for the second order formalism, the first order formalism version can be obtained in a similar way. Firstly, 
let us define the variables; we have two equations: the decomposition formula and the Gauss-Codazzi equation 
written compactly in (13). Since we are working in (2-|-l)-dimension, it is natural to use matrix group SO{3) to 
represents the 3-dimensional intrinsic curvature, and the subgroup SO{2) to represents the 2-dimensional intrinsic 
curvature. But for simplicity of the calculation, we use the unitary group SU{2) -the double cover of SO{3)- instead 
of SO{3). Having information of G SU{2) is equivalent with having information of the plane and angle of the 3- 
dimensional rotation, which are denoted, respectively, by {J G su( 2 ), dd G K}. The same way goes with G SO(2), 
it contains same amount of information with the plane-angle of the 2-dimensional rotation {j G 20 ( 2 ), Se) G K}. Using 
the exponential map, we can write: 

= exp {jS(j)) = 72x2 cos S(j) -G j sindcj), (14) 


with 

0 r 

[ -1 0 J ’ 

is the generator of SO{2). Doing the same way to the element of SU{2), we obtain: 

SO so 

~ exp {J59) = 72x2 cos -I- Jsin 

with J G su( 2 ) is an element of the Lie algebra of SU(2), satisfying: 


J = Ada, 


rJa = 1 . 


o-Q are the basis of su( 2 ), namely, the generator of SU{2): 


(15) 


(16) 


satisfying the algebra structure relation as follow: 




( 17 ) 
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FIG. 1. Suppose we have a 3-dimensional curved manifold Ma discretized by four tetrahedron in the figure above (in flat case 
it is known as 1-4 Pachner move). Then we take an embedded slice Ea as the surface of one tetrahedron (the dark blue surface 
discretized by three triangle). Embedded on Ea, we take a loop 7 circling a point of a tetrahedron. Attached to 7 , are the 
SU{2) and 50(2) holonomy, which are related to the 3D and 2D intrinsic curvature, respectively. 


with (To are the Pauli matrices: 



■ 0 

1 ■ 


' Q -i' 


r 1 0 1 

(J^ — 

1 

1 

0 

, CTy — 

i 0 

II 

b 

— 1 

t-H 

1 

0 

_1 


The ^ factor in (16) comes out from the normalization occuring when we write SO{3) using SU{2) representation, 
i.e., to reduce the factor 2 in relation (17) 

Taking the trace of (13) gives the relation between 3-dimensional and 2-dimensional rotation angle: 

tr^iJ = tr^Tf -f trT, (18) 

2 cos ^ 2 cos 54 > 


using the fact that the elements of the algebra are skew-symmetric. The other relation we need to have the full 
information contained by (13) is the planes of rotation. Since the embedded surface is 2-dimensional, the plane of 
rotation for Scj) is trivial, which is (15). The plane of rotation J for 69 can be obtained by solving complex 2x2 
matrix linear equation (16): 


J = 


^H-L 


se 


2 x 2 cos 


sinf 


(19) 


In conclusion, the contracted Gauss-Codazzi equation concerning the relation of curvatures of a manifold and its 
submanifold in (2-1-1) dimension can be written using holonomy representation (13) or using the plane-angle repre¬ 
sentation, i.e., the relation between rotation angles by (18) and the condition for the planes of rotation, by (15) and 
(19). 


III. DISCRETE (2+1) GEOMETRY 
A. Geometrical setting 

Using the Gauss-Codazzi equation in terms of holonomy and plane-angle representation obtained in the previous 
section, we can write the Gauss-Codazzi equation for discrete (2+1) geometry; this is the objective of this section. 
Let a portion of curved 3-dimensional manifold M be discretized by flat tetrahedra, we call this discretized manifold 
Ma. See FIG. 1. 


1. Angle relation 

In this section we derive the Gauss-Codazzi equation purely from the angles between simplices, i.e, the relation 
between geometrical objects, without reference to coordinates (this is, indeed, the reason Tullio Regge developed his 
Regge calculus and discrete general relativity, in his famous paper titled ’General relativity without coordinates’ [5]). 
There are two kinds of relevant angles in a flat tetrahedron: the ’angle at a vertex’ between two segments of the 
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FIG. 2. Given angles tjiab, (j>bc, <j>ca at point p of a tetrahedron, we could obtain the dihedral angle Oabac- In fact, Oabac is only 
<j>bc projected on the plane normal to segment a. 



FIG. 3. Suppose Ea is a 2D surface discretized by three triangles. Point p is the hinge where the deficit angle Sep is located. 
S(f) is the 2D intrinsic curvature of Ea at point p. 


tetrahedron meeting at a vertex, which we denote (j)] and the ’dihedral angle at a segment namely the angle between 
two triangles, which we denote 9. These two types of angles are related by the ’angle formula of the tetrahedron’: 


cos Ofibac 


COS (pbc - COS 4>ab COS (pac 
sin (pab sin epac 


( 20 ) 


See FIG. 2. 


2. Isosceles tetrahedron 


For simplicity of the derivation in the next section, we consider isosceles tetrahedron. An isosceles tetrahedron 
is build by three isosceles triangles and one equilateral triangle. In this case, all angles around one point p of the 
tetrahedron are equal, say (pi = (p, and the dihedral angle relation (20) becomes: 


cos 6 = 


cosep 

(1 + COS(p) ' 


This implies that the dihedral angles between two isosceles triangles are equal as well, say 9i = 9. Another geometrical 
property of a tetrahedron which is useful for the derivation in this section is the volume. The volume of unit isosceles 
parallelepiped is: 


vol = (ni X n 2 ) ■ ns 


i (1 — cos(p) \/l + 2 cos (p 



( 21 ) 


B. Curvatures in discrete geometry 

In n-dimensional discrete geometry, the notion of curvature is represented using the plane-angle representation. 
The intrinsic curvature is represented by the angle of rotation, which in this case is called as deficit angle, and this 
deficit angle is located on the hinge, the (n — 2) form duals to the plane of rotation. The geometrical interpretation 
of deficit angle is illustrated by the FIG. 3. 

To obtain the 2-dimensional and 3-dimensional intrinsic curvatures in discrete geometry picture, firstly we need to 
define a discretized surface embedded in Ma, say Ea. The simplest way is to take the surface of the tetrahedron as a 








p 



FIG. 4. A portion of a discretized curved manifold Ma with an embedded discretized hypersurface Ea- It is clear that the 
loop 7 circles point p, where the 2D deficit angle 54) is located. S4> represents 2D intrinsic curvature since it is located on 
a point. The loop also crosses 3 segments which are si, S 2 , S 3 . Each segments contains its own 3D deficit angle 56, which 
represent 3D intrinsic curvature since it is located on a line. Using the exponential map, we can obtain holonomies on each 
line corresponding to the deficit angles. From these holonomies, we can obtain the 3D total holonomy along loop 7 , which is 
the product of holonomies on the 3 segments. 


slice (see FIG. I). Having Ea embedded on Ma, we consider a loop 7 on Sa. We consider the holonomy related to the 
2D and 3D curvatures defined by this loop. The 3D curvature is defined by the 3D holonomy S 5'0(3) ^ SU{2) 
around loop 7 , which describes the intrinsic curvature of Ma, while the 2D curvature is the 2D holonomy G 50(2), 
describing the intrinsic curvature of slice Sa, still along the same loop. Taking trace of the holonomies, we obtain 
the angles of rotation, i.e., the amount of a components of a vector parallel to the plane of rotation get rotated when 
parallel transported along loop 7 . In discrete geometry, this rotation angle is the same as the dehcit angle on the 
hinges. 


1. 2D intrinsic curvature (SO(2)) 

Using FIG. 4 as the simplest case, we take loop 7 circling point p on the surface E, therefore the 2D curvature is 
the deficit angle on p: 


^R = 6(j) = 2TT-'^, 


or in terms of the cosine function: 


cos 


3 3 

S(j) = cos (j)i — cos (pi sin pj sin • 


2 = 1 


For the isosceles tetrahedron case, (22)-(23) become: 

^R = 5(I) = 2tt- 3(j), 


( 22 ) 


(23) 


cos S(j) = cos^ (/) — 3 cos (j) sin^ (f = cos if (4 cos^ if — 3) . 

We could obtain the holonomy by the exponential map (14), using the 50(2) algebra (15). 


2. 3D intrinsic curvature (51/(2)) 

We could write 3D curvature by the holonomy of SU(2) using (16), but the 3-dimensional holonomy along loop 7 
circling point p is a product of three holonomies on the three hinges, because the loop crosses three hinges. See FIG. 

4. 

2=1 


( 24 ) 
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Using (16), we obtain: 

3 

r I T rns ^ 


TT / T T ■ 

77 = 11 ( / COS — + Sin — 


i=l 


2=1 


COS 


S9^ 

f ^ J.sin^cos Ycos^ + 

i—l,i^j,k,j<.k 


J . S0i . 50j 50k T-r 

JiSin—JjSin^cos— +11 Ji 

i—l,i<.jyi,j^k i—1 


Sin 


50^ 
2 ’ 


where each hinge (segment) i is dual to plane Ji = Ji ■ a and is attached by a deficit angle 50i. Taking trace of the 
equation above gives: 


tP^ =Z^llcos—+ 2 _^ trdi sin — coscos—h tr (J^ Jj) sin — sincos ^- 

2cos^® 2 *=1 i=i,i^j,k,j<k q i=l,i<j,i,j^k 


- 2 (J,.+ ) 


tr 






sm - 


^2 = 1 / 2=1 


or simply: 


+ 2 ^ Ji X J 2 )■'/a 


50 A 5^, 

cos—= ||cos—- 
2=1 


y~! (ji • 


2 = 1 , 2 <J, 2 ,J 7 ^/C 


— sin —^ 
2 2 


sin — sin —^ cos-h 


50 k 

2 


(/i X J 2 ) 


J 3 sin 




But • Jj^ = cos^i, and ^Ji x is just the volume of the parallelepiped defined by | Ji, J 2 , (Remember 

that J = *J, or 1-form J is the Hodge dual of the 2-form J in 3-dimensional vector space. The direction of J is 
normal to the plane discribed by J). Therefore, the 3D intrinsic curvature is:^ 


50 -p-r 50i i • 50i . 50j 50k , ,ntt . 50i 

cos—=||cos—- Z^ cos((i,sm —sm^cos—+ (vol)||sm—, 

2 = 1 2 =l, 2 <_ 7 , 2 ,J' 7 ^fc 2=1 


(25) 


or 


3 r An o I TT X • 50fc , ,, tt . 50^ 1 

7t = d0 = 2 arccos I I cos —- > cos (pi sm — sm -y cos ——[- (vol) I I sin — , 


I Z=1 


i—l,i<.j,i,j^k 


2=1 


with vol is the volume form described in (21). For isosceles tetrahedron case, we have (pi = (p and 50i = 50, so the 
holonomy is: 


50 


50 . 50 


T 3 c (7 / o 0(7 . oe \ V 

77 = 7 cos^ y + cos y sm — ) ^ J, 


50 . ,50 

O' + I cos — sm — 


vi=l / '' ' ij = l,i^j 

Therefore, the 3D intrinsic curvature for isosceles tetrahedra case is: 


Y, • < 7 ) (Z • + sin^ y n > 


3 I 3 50 50 . 2 50 / ,N . 3 50\ 

77 = 00 = 2 arccos ( cos — 3 cos (p cos y sm y + (vol) sm y ) • 


5. .8D extrinsic curvature 

Equation (18) is a linear relation between the trace of 3D holonomy with the 2D holonomy and its ’residual’ 
part T. Inserting (23) and (25) to (18), we could obtain tiT, but this is not the extrinsic curvature, since T = 


^ (25) has a similar form with (23), except the cos pi anci the term which contains the volume form. 
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FIG. 5. From the dihedral angle relation, we could obtain the internal dihedral angle 6^ on a segment, a is the external dihedral 
angle, while 9' is the angle between normals of the two triangles. 


S + [K, K] + {I 3 — I 2 ) ■ Next, we would like to obtain an angle describing the extrinsic curvature K. We do not expect 
the relation between this angle with and to be linear, since linearity is held on the holonomy representation of 
Gauss-Codazzi (13), while the relation between their rotation angle ^R and ^R does not need to be linear. 

Consider the isosceles tetrahedron case. On each segment of the tetrahedron, lies the 3D deficit angle 50. This 
deficit angle is defined as: 


(50 = 27r — (0 + a), 

where 9 is the internal dihedral angle, and a is the external dihedral angle coming from the dihedral angles of other 
tetrahedra (remember that the full discretized (2 + 1) picture is the 1-4 Pachner moves). Let’s introduce the quantity 

9 = a — TT. (26) 

For the case where is flat, 59 = 0. This causes 0 -I- o; = 27r, and using definition (26), we obtain: 

9 + 9' = Tr. 

In this flat case, it is clear that 0' is the angle between the normals of two triangles, see FIG. 5. 

Therefore, 9' is in accordance with the definition of extrinsic curvature in (3), where K is defined as the covariant 
derivative of the normal n to the hypersurface S. Because of this reason, we define 0 as the 2D extrinsic curvature, 
since in a general curved case, it will inherit the curvature of the 3D manifold. 


C. Discrete Gauss-Codazzi equation in (2-|-l) dimension 

A simpler but equivalent way to obtain the 3D instrinsic curvature is to use the purely geometrical picture of the 
discretized manifold. See FIG. 6. 

From Figure 6, we obtain the projected 3D intrinsic curvature of each segment on the ’artificial’ segment dual to 
the total plane of rotation as: 


Using 


2 <50L sin^x'cos^^ 

2 1 — COS^ x' COS^ ^ 

Ji = siny'fc + cosyfi. 


(27) 


we obtain: 


(^Ji, Jj'^ 


= cos 


sin^ x' + cos^ x' cos^ 


27r 

Y’ 


2 

cosY' = 3 (1 “ ’ 

sin^ y' = ^ (1 -I- 2 cos (j >), 
o 


or: 


(28) 

(29) 
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FIG. 6. The 3D deficit angle on each segment can be thought as a ’piece of space’ inserted along the hinge. The ’size’ of this 
’piece of space’ is described by the dehcit angle S6, which is located on the plane normal to the hinge. The total deficit angle 
50 obtained by adding three rotation in (24), is located in the plane normal to an ’artificial’ hinge in the direction k, since we 
use isosceles tetrahedron. Therefore, projecting SO to this plane, we obtain 59\^. 


since we use isosceles tetrahedron. All the curvatures: 3D intrinsic, 2D intrinsic, and 2D extrinsic curvatures are 
contained in (27), through: 




T 




(30) 


27r — 5(j) 27r — 


(31) 


59 = TT — {9 + O') = TT — 


^arccos 


coscj) 

1 + cos (/) 



(32) 


Using the half-angle formula and inserting (28)-(32) to (27), we obtain the relation between angles of the discrete 
Gauss-Codazzi equation for a special discretization using isosceles tetrahedra: 


1 — cos K -b arccos 


cos 


2-n-R 


1 + cos 


1 -b cos 


1 - I (1-cos^) (1- 


cos 


(33) 


Together with their planes of rotation: j in (15) for ^R and k ■ a for (since we use isosceles tetrahedron), (33) is 
the Gauss-Codazzi equation in plane-angle representation, which contains same amount of information provided by 
the linear Gauss-Codazzi equation in holonomy representation in (13). 


IV. SPIN NETWORK STATES AND 3D LOOP QUANTUM GRAVITY 

A. 2-complex and its slice 

In this section, we apply the results obtained above to loop quantum gravity (LQG). We will briefly review the 
quantization of gravity in first order formalism using the loop representation. From the real space representation, 
we move to the (Hodge) dual space [7]. (The dual space is introduced in analogy to what is usually done in higher 
dimension.) In LQG, the fundamental geometrical object is the 2-complex. A 2-complex is a 2-dimensional geometrical 
object dual to the discrete manifold Ma. Let us call the 2-complex as *Ma. We can take a slice on a 2-complex, which 
is a 1-complex called as *Sa, dual to the discrete hypersurface Ea of A 1-complex is a graph consisting oriented 
lines called links I and points called nodes n. The variables attached to the graphs are the variables coming from 
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FIG. 7. The figure above is a 2-complex dual to the geometric figure in FIG. 1. Suppose we have a 3-dimensional curved 
manifold Ma discretized by four tetrahedron as the first figure above (in flat case it is known as 1-4 Pachner move). Then we 
take an embedded slice Ea as the surface of one tetrahedron (the dark blue surface discretized by three triangle). Going to 
the dual space picture, Ma is dual to a 2-complex *Ma called as the 3D bubble [7] (its geometry is defined by black lines and 
black points), while Sa is dual to a 1-complex *Ea called as 2D bubble (defined by purple lines and point). Glearly, the 2D 
bubble *Ea is embedded on the 3D bubble *Ma. The bubble graph define complete curvature of the portion of space. 


first order formalism: the holonomy Ui G SU(2) and algebra Ji G su( 2 )* on each links I of the graph. A graph with 
elements of SU{2) x su( 2 )* is called a spin network. In the canonical quantization, we promote {Ui, Ji} to operators 

|l7i, Ji|. The space SU{2) x su( 2 )* is the phase-space of LQG, and the pair Ji| are the pair conjugate to each 
other, satisfying a quantized-algebra structure as follow [3, 7, 11]: 


Ui,Uv 


= 0 


=-6u'a^Ui 


fi p 


= -SwelJt. 


On each node, the gauge invariant condition must be satisfied: 

^ 1^) = E 1^) = 0- 


iGn 


This condition selects an invariant Hilbert space JC = L 2 [SU{2y/SU{2)'^] 9 
I A, A, A| [3, 7, 11]. As an example for our case, see FIG. 7. 


of a closed triangle formed by 


B. Curvature operators 


Given a 2D bubble graph with a phase space variable on each link, we can construct the operators of 3D, 2D 
intrinsics and 2D extrinsic curvature from these variables. 

We construct these operators using the basic phase-space operators, for example, the angle operator from the 
phase-space variables [ 12 ]: 


(pi = arccos 


{ A 

2 

+ 

Jk 

2 

A 

2 

V ^ 

A 

Jk 



i, j, k Gn. 


From this angle operator, we obtain the 2D intrinsic curvature along loop 7 as: 


2i? = 27r-EA, 

2 C7 
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{Us^Js} +ZA 

a^,j 2 } 



FIG. 8. Taking only the slice on the dual space representation, we obtain the spin network *Ea. 7 is the loop chosen. Inside 
7 is where the point p located in the real space representation. 

which is clearly the deficit angle on a point in the real space representation. See FIG. 8. 

The 3D intrinsic curvature can be obtained from the holonomy around a loop. Remember that the direct geometrical 
interpretation of spacetime can only be obtained from second order formalism, while in LQG, the fundamental variables 
comes from first order formalism: the holonomy U comes from the curvature 2-form F, instead of the Riemannian 
curvature R (see Section 1). The relation between F and R, subjected to the torsionless condition is: 

R = e{F), (34) 

with e is the local trivialization map between them. Gontracting (34) with an infinitesimal area inside the loop (carried 
by the derivation in Section 1) gives: 

H = e{U), U = e-^{H), 
since e is a diffeomorphism and must have inverse. 

The 3D intrinsic curvature is ^R ^ tri/, with F[ is the holonomy coming from the contraction of Riemanian 
curvature R with an inhnitesimal area. But since e is a diffeomorphism and trace of the holonomy is invariant under 
diffeomorphism and gauge transformation, we obtain triJ = tvU. Therefore, we can write the 3D intrinsic curvature 
operator as: 


^R = 2arccos 


trC/iC/2t/3 


(35) 


The extrinsic curvature operator can be obtained directly from the discrete Gauss-Codazzi equation (33): 


K = arccos 1 — 


1 + cos 


cos 


2 ir-^R 


1 - i (l - COS ( 


1 - COS 


— arccos 


1 + cos 


C. The semi-classical limit 


In this section, we show that the linearity of the relation between tc^Fl and tr^H in equation (18) is recovered 
from the spin network calculation in the semi-classical limit, by dropping the Maslov phase. Let us take a simple spin 
network illustrated in FIG. 9. For simplicity in our calculation, we use a closed graph as an example, but the result is 
valid for any spin network graph. The vector and covector state of the spin network *Sa in the algebra representation 
are: 


1 ^) = Y, \j^rnini ... jemeUe), 

V ni..nQ ) 


(36) 


E 


•rnf-rrir,n<-, •rne^rrir.ri-x n 


3i<« ... j^TOgngl , 


fi/ / 


(37) 


with i is the intertwinner attached on the node. See [7, 11, 13] for a detail explanation about the spin network state. 
Ghanging the representation basis using: 


{U\jmn) = {U) 
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FIG. 9. A simple spin network *Sa (the graph in pnrple), which will give a geometrical interpretation of a ’fuzzy’ quantum 
tetrahedron in the semi-classical limit. This spin network can be obtained from the spin network in FIG. 8, by identifying 
their open links together. In 3D theory, the nodes (purple points) are dual to triangles, and the links (purple lines) are dual to 
segments. 


with (U) is the matrix representation of C/ S SU(2) in (2j + l)-dimension, we can write the vector and covector 
states in the group representation basis: 


= E 


dUi..dUe 


i=l 


r ® 

(l/;| = ^ -n'n'n' / . .df/g „,([/') (C/( . .t/g | . 

(^ J i=l ' ' 

\ n^/ ..Tig/ / 

We choose a specific closed path on spin network *Sa) say, loop 7 , defined by links {1 — 2 — 3} (See FIG. 9). The 
holonomy U = U 1 U 2 U 3 are attached along 7 , therefore we have a well-defined operator related to the 3D intrinsic 
curvature operator (35) on the loop: 


tT^H= =tT{UiU2U3)=J2SmnDl„{UlU2U3). 


2 cos 


Since the matrix representation satisfies: 

(j, m\ D (C/iG^t/a) |j, n) = (GiC/aC/a) = E (7) {U 2 ) Dl (Us ), 


kA 


we can write 


tr^H= Y. 5^nDl^{U,)Di,{U2)Dl{U,). 


(38) 


(39) 


(40) 


k.l', 


The geometrical interpretation of the spin-} in (38)-(40) is the measure on the ’artificial’ hinge dual to the plane of 
rotation inside loop 7 . See FIG. 6 and Section III G. 

Acting tr^H with the spin network states gives: 


\tl= I j^mem3n2j^mem'^n2j^m3lrn'^^n2ln'2 j j 

vmQ,l,m^,n 2 ,m^/,n 2 / j \ ms,fc,m2,ni, 7712 /,' 11 . 5 / 


^m5m2ni^m5m2ni^m2A:m2^niA:ni . ^ 


... X 


( 'y jrnirnin3^mirn'in3jn3Tnn'3jrnirnm[ j j y ^ 

m 4 ,m,mi , 713 ,mj^/, 7 X 3 / / \ 7 i 4 .. 7 i 6 




but the first three-terms in the parantheses are only the Wigner 6 j-symbols: 


' ^ '' Js 72 J / I 72 7 } 7/1 j'l 73 7 




(41) 


nLiisii. 
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FIG. 10. 
which is 


The 6 j-symbol 


can be represented using a graph (left), where it can be translated 


J6 J3 J2 

is is i 

a tetrahedron. Node a on the graph is the triangle a on the tetrahedron. 6i is the dihedral angle 


into its dual (right) 
on j. 


We are ready to take the semi-classical limit by setting all the jiS to be large. Using the result proved by Roberts 
[14, 15], for large spin j, the 6 j-symbol can be approximated as follow: 


{ 6 i} = 


Jl J2 J3 

ji is ie 




cos 



where the six constructs a tetrahedron, Oi are the internal dihedral angle on segment ji, and jUj is the volume of 
the tetrahedron. The j term is known as the Maslov phase. Using this result to our case, we obtain: 

« + 2i303 + 2i26»2 -b i 6 »/ - 0 . (42) 

where the last step we use the fact that ji = j'i, noting that j' comes from the covector state (37). This fact also 
cause 9i = 9'i. See FIG. 10. Remember that j is the measure on the ’artificial’ hinge inside loop 7 , in other words, 
the spin number j is the ’length’ of the artificial hinge. Along this artificial hinge, there is a point p where the 2D 
deficit angle is located. Taking ji ^ j, we have: 


9q « 0, 03 « -, 02 ^ 9i ^ (j)23 — </’6) 


See FIG. 10. Using this approximation, (42) becomes: 

1 


J6 J3 32 

‘J A' 


■^3 ^2 j J v^l27r|Ui (j6,j3,j2, j)| 

Doing the same way for the other two 6 j-symbols, we obtain: 


35 32 3i 

'J A' 


■^2 Jl J J \/l27r|U2 (js, j2,Ji, j)| 


cos (^(j3 + J 2 ) TT + j^5 
cos ((j2 + 71 ) 71 +j(('5 



COS (?1 

\/l 27 r|U 3 (74,71,73, j)| V 

By writing the cosine function using the Euler formula: 

cos ^ = 2 ~*a;) 

we can write the products of the three 6 j-symbols as: 


f 74 jl 73 

1 fi 73 3 


+ 73)7r+7()i4 



n“ 1 i2„Vi2>IUIIUIIC,| " (™ +*+*)-+7- t) + - 

... + cos ^ 2737 r + j (^4 — 4>5 + (je) ~ 4 ) (^7271 + j {—(j)4 + ((>5 + ipe) ~ ”*"■■■ 

... + cos (^ 27 i 7 r -bj {(j )4 + <^5 - <^6) - 0 ) ■ 


(43) 
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From this point, we take the sum of ji +j 2 +j 3 to be an integer/natural number (we discard the half-integer possibility) 
and choose the spin number j = 1, then we could write: 


n 


2=1 


1 1 


cos ( S(j) 



(44) 


where we denoted the remaining terms by O and using (j )4 + + (j)Q = 2Tr — 6<j) (see FIG. 9). Inserting (44) to (41), 

we obtain: 



The linearity relation (18) can be recovered by dropping the Maslov phase: 

An important point is we have a freedom to choose the spin number j, since this spin number comes from (38), 
which is the dimension of the SU{2) representation. If we set j = 0, then, the volume of the 6j tetrahedra \ Vi\ will be 
zero and relation (43) will diverge, given any j's. But since the range of spin j is j > with j multiple of half, then 
the spin number j acts as an ultraviolet cut-off which prevents relation (43) from divergence. 


V. CONCLUSION 

We have obtained the Gauss-Codazzi equation for a discrete (2-|-l)-dimensional manifold, discretized by isosceles 
tetrahedra. We have studied the definition of extrinsic curvature in the discretized context. With definitions of the 
3-dimensional intrinsic, 2-dimensional intrinsic, and 2-dimensional extrinsic curvature in the discrete picture, we have 
promoted them to curvature operators, acting on the spin network states of (2-|-l)-dimensional loop quantum gravity. 
In the semi classical limit, we have shown that the linearity between tr^H and tr^iJ in equation (18) is recovered in 
the spin network calculation, by dropping the Maslov phase. There exist a natural ultraviolet cut-off which prevents 
the discrete Gauss-Codazzi equation from divergences in the semi-classical limit. 
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